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Abstract 

Spectrum of the Laplacian on spherical domains is analyzed from the point of view 
of the heat equation on the cone. The series solution to the heat equation on the cone 
is known to lead to a study of the Laplacian eigenvalue problem on domains on the 
sphere in higher dimensions. It is found that the solution leads naturally to a spectral 
function, a 'generating function' for the eigenvalues and multiplicities of the Laplacian, 
expressible in closed form for certain domains on the sphere. Analytical properties of 
the spectral function suggest a simple scaling procedure for estimating the eigenvalues. 
Comparison of the first eigenvalue estimate with the available theoretical and numerical 
results for some specific domains shows remarkable agreement. 

The heat equation on the n-dimensional cone has been addressed by various authors 
in the past and series solutions have been obtained. The n = 2 solution was obtained 
by Sommerfeld [1894]. The n = 3 case was considered within the context of circular 
cones by Carslaw and Jaeger [1959]. For higher dimensions, the applicable solution has 
been presented by Cheeger [1983]. The radial component of the series solution involves 
modified Bessel functions while the angular component involves the eigenvalues and the 
eigenfunctions of the Laplacian on a domain on the n—1 dimensional sphere. It is found that 
the solution leads naturally to a spectral function, a 'generating function' for the eigenvalues 
and their multiplicities, expressible in closed form for certain domains on the sphere such 
as the octant triangle on the two-sphere and analogous ones on higher dimensional spheres. 
Analytical properties of the spectral function suggest a simple scaling procedure to estimate 
the first few eigenvalues for related domains. The estimates for some specific domains are 
found to be in excellent agreement with the available theoretical and numerical results. 



1 Heat Equation on the Cone 

Consider M" with coordinate vector x, conveniently split into radial and angular parts, r 
and f (when r 7^ 0), as 



^2 



Unit radial vectors f trace out the n—1 dimensional sphere /S""^ at = 1. Consider a 
domain O""^ of dimension n — 1 on S""^ with boundary 5^7"^^ (that can be empty). Let 
Z?" be the (unbounded-) cone through 17"^^: the region of M" comprising of all vectors x 
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whose unit radial vectors f are in ^ (for definiteness, r = is included). Its boundary 
dD^ is similarly the cone through d^^~^. Now, consider the heat equation on the cone D", 

i ' 

is the Laplacian, r is the time variable and f{x,x',T) is the heat kernel required to 
satisfy the initial condition f{x,x',0) = 6{x — x') and a boundary condition on dD^. 
Dirichlet boundary condition requires / to vanish on dD^. Neumann boundary condition 
requires the normal derivative of / to vanish on Standard arguments involving the 

Laplacian show that the solution to the heat equation is unique under these conditions and 
f{x,x',T) regular at r = going to zero sufficiently fast as r — 5- oo. 

An important property of the heat kernel is that it is multiplicative in a Cartesian 
product. Consider a cone D"^ that is a Cartesian product D^^ x of two cones D^^ and 
where n = ni + n2, formed by pairing the coordinate vectors as a; = (xi, a;2). Then it 
follows from the heat equation that the heat kernel / on is given by 

f{x,x',T) = fi{xi,x[,T)f2{x2,X2,T), (3) 

where fi and /2 are the heat kernels on D^^ and respectively. The initial condition 
at r = is satisfied. As can be verified, Dirichlet or Neumann boundary conditions on 
D^^ and imply analogous or mixed ones on (that consists of dD^^ x and 

X dD^^ joined along dD"^' x dD'^'). 
There are some special cases when the heat kernel can be explicitly written down. When 
the cone is M, that is the whole real line, the solution is 

f{x,x',T) = ^e-M--^'^\ (4) 
v27rT 

that tends to 5{x ~ x') when r — t- as required. When the cone is M>o, that is the 
nonnegative real line, with the origin x = considered as its boundary, the solution is 
obtained by the method of images. 



/(x,x',r) = [e-M^-^'f - e-M^-^^')') = \f^erM^'^^'') sinh (^) . 



(5) 



This is for Dirichlet boundary condition. For Neumann boundary condition, we will have 
cosh in place of sinh above. 

These solutions can be combined to obtain solutions on their Cartesian product spaces 
making use of the multiplicative property of the heat kernel mentioned above. We thus 
have the heat kernel on as the product of n copies of (4), 

/(x,x',r) = -^e-^(="--')'. (6) 
(27rr) 2 

We also have the heat kernel on R>q as the product of n copies of (5) (under Dirichlet 
boundary condition), 

fix,x',r) = (^)%--(^^^^'^)n--l^(^) • (7) 
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We can also combine the solutions, say p copies of (4) and q copies of (5), to obtain heat 
kernels on product spaces such as MP x M>q. 

However, the heat equation on the cone is not solvable in general, except as a series 
solution. One of the ways of approaching the series solution is to represent it as 

/(cc, cc', r) = V / df^tJ,e~^^'^''(j)^,ua{x)(j)^ua{x'). (8) 

ua 

Here (pnua^s are the eigenfunctions of the Laplacian on the cone of eigenvalue — 
having additional labels ua for later use, 

V"^ (Pij.ua = -IJ-'^'Pfjua- (9) 

It can be verified that representation (8) does solve the heat equation, (p^u^s are expected 
to behave similarly to the heat kernel on dD^, at r = and as r ^ oo. Standard arguments 
show that we then have /x^ > and that ji^ covers all of the positive real line because of 
the unboundedness of the cone. Initial condition will be satisfied with the (p/xua^s forming a 
complete system, normalized according to 

/ (r-X(p^t,a{x)(pfj.'u'a'{x) = -d{ll - l^')5i,iy>6cra'- (10) 

We thus have the Laplacian boundary value problem (9) to be solved that can be rewritten 
in radial and angular coordinates involving the spherical Laplacian. 



2 Spherical Laplacian 

Going to radial and angular coordinates, one obtains the following expression for the Lapla- 
cian on the cone acting on some function (p, 

^^*-^^^(^"-'f)->^* 

where V|. acts on the angular dependence. Identifying the angular coordinates as the 
coordinates on the sphere S*""^, one notes that V|. is the Laplace-Beltrami operator, or 
simply the Laplacian, on the spherical domain Q""^ of interest. 

Angular dependence can be expressed as a dependence on r = ^ so that functions on 
can be viewed as zero-degree (positive-)homogeneous functions of x in the cone D"' 
(a scaling of x by say a scales r by \a\ so that functions of r remain unaffected if a > 0). 
V5 can hence be identified with r^V^ acting on zero-degree homogeneous functions of x. 
Of interest to us is the Laplacian eigenvalue problem on J7"~"^, 

VlKaif) = -XK^{f), A = i/(i/ + n-2), (12) 

where the zero-degree homogeneous function hi^^{r), a labeling any multiplicity, is the 
eigenfunction satisfying appropriate boundary condition on dil'^~^. As can be verified, this 
eigenvalue equation can be equivalently expressed as the Laplace equation V^(r'^/ii,o-) = 
in i)" for a i^-degree homogeneous function r'^hvair)- It is thus natural to work with the 
degree variable u in place of the eigenvalue A. 

Boundary value problems of the above kind have been extensively studied and it turns 
out that the eigenvalues are, and hence the z^'s can be taken to be, real, non-negative. 
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unbounded and discrete, and that the eigenfunctions can be taken to be real and form 
a complete system. We will assume that the eigenfunctions are normalized to form an 
orthonormal system 

[ (r-HK^{f)K,^,{r) = 5,,>5^„,, (13) 
where (P~^r denotes the volume element on the unit sphere S"'~^ (area element if 5^). 



3 Series Solution on the Cone 

Going to radial and angular coordinates involving the spherical Laplacian as in (11), the 
Laplacian boundary value problem (9) on the cone can be rewritten as 

where V| is the Laplacian on $7""-^ discussed above. The complete set of eigenfunctions 
hua{r) of V|, of eigenvalues + n — 2) and multiplicity label a, enables us to look for 
4'^.i'a{x) of the form 

n — 2 

(t>fiua{x) = r ~g^ua{r)Ka{r). (15) 

Boundary condition for ^fiu<j{x) on implies analogous one for hi,a{r) on Ofi""^. Nor- 
malization (10) can be ensured with the g^ua'^ normalized according to 

noo ^ 

/ drrg^^a{r)g^'va{r) = -5{iJi-ij!). (16) 
Jo A* 

The action of V| on hy(^{f ) is given by (12) so that (14) gives rise to 



This is the Bessel's differential equation with the acceptable solution g^vair) = J^_,_2i=2 (/xr). 
Behavior of gfiuaif) as r — )■ and r — )■ oo is consistent with the expected one for (j)^i,u{x). 
Orthonormality and completeness of g^uaif) follows from the properties of the Bessel func- 
tion as the basis of the Hankcl transform. We thus have 

n — 2 

4>iJ.va{x) = r 2" J^+rt=2 (iir) Ka{r)- (18) 
Using this in the heat kernel representation (8), we get 

2 1 

U 2 2 ^ 

Making use of an identity that expresses the integral over product of two Bessel functions 
in terms of the modified Bessel function, we obtain the series solution 



f{x,x',r) = ^{rrr'^e-Mr'+^"^Y.Iu+^ (^) /i..(r)/i..(r'). 



(20) 



This result was obtained under different contexts by various authors. For n = 2 it was 
obtained by Sommerfeld [1894]. For n = 3, it was considered within the context of circular 
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cones by Carslaw and Jaeger [1959]. For general dimensions, it has been presented by 
Cheeger [1983]. The leading term in the series (20) can be obtained by making use of the 
expansion for the Bessel functions, 

f{x,x',T) ~ I f!fY e-Mr'+r'')h.^{r)K^{f'), (21) 

r (z/1 + f) (2t)2 \2tJ 

where i^i is the first ly and F is the Gamma function. 

In the special case of = 2 in polar coordinates r and 9, we have the cone < ^ < </? for 
some <2'K having the circular arc at r = 1 as the spherical domain. For the eigenfunctions 
huai'f^) under Dirichlet boundary condition, we have ■^J^sm.{v6) with multiplicity one and 

= = 1, 2, • • • (under Neumann boundary condition, we have instead .^/^ cos(y9) and 

also for A; = 0). For </? = tt, i^'s are positive integers, and for <^ = |, they are positive 
even integers. In these cases, the heat kernel is known explicitly and the series solution is 
consistent with the known identities involving modified Bessel functions. 



4 Product Domains 

As noted in section 1, an important property of the heat kernel is that it is multiplicative in a 
Cartesian product. Because the series solution to the heat equation involves the spectrum of 
the spherical Laplacian, the multiplicative property has implications on spherical domains. 
Towards this end, it is helpful to have a look into product domains on the sphere implied 
by the Cartesian products. These product domains are recognizable as topological 'join's 
of the participating domains. In our case, they can be constructed as follows. 

Let on S*"!"^ and ^2^""^ on S'^'^^'^ be two spherical domains. Consider their cones 

L>"i in R"i and D'^'' in . Let the product cone D'^^ x Dl^' in M"i+"2 be the Cartesian 
product of I?"^ and D^"^ . The spherical domain associated with D^^ x Z)" ^ then gives us 
the product domain *Q2^~^ on 5'"i+"2-i obtained by restricting the combined radial 

coordinate rf + r\ to unity, il"^""^ * ^^^^^ nontrivially related to and Vi^^^^ ■ Its 

dimension ni + n2 — 1 is one more than the sum of the dimensions of the two participating 
domains. Its boundary consists of 9^7"^^^ * and fi"^"^ * d^^^"^ (one or both of which 

could be absent) joined along their boundary * d^^^^ . 

Let denote the 'zero-dimensional sphere': the two-point set {-|-1,— 1} having its 
'cone' as M, the whole real line. Also, let be the one-point set {-1-1} having its cone 
as M>o, the nonnegative real line. Then, as can be verified, * abbreviated as 5°*^ 
is simply the unit circle . More generally, S"""^ can be obtained as 5"*^*" by taking the 
♦-product of n copies of . Similarly, it can be verified that = T^*"^ is the quadrant arc: 
the 90°-arc of the unit circle occupying one-fourth of the circumference. = 2^0*3 jg ^.j^g 
octant triangle: the spherical triangle on the two-sphere having three 90° angles occupying 
one-eighth of the spherical surface. = T^** is the spherical tetrahedron on occupying 
one-sixteenth of the sphere. More generally, T"^^ = T''*" is a spherical polytope on S*""^ 
occupying (2"'*)*'^ of the sphere. We can also construct *— products comprising of both 
and T°. For instance, 5° * T° is the half-circle, * T° is the half-sphere or half-S"^, and in 
general S"'~'^*T^ is half-S""^^. More generally, SP~^*T'i~^ is a domain on 5^+9-1 obtained 
by cutting away 5*^+^^^ perpendicularly into half g-times, for instance following the steps 

^p+g-l _^ gp+q-2 ^ ^ gp+q-i ^ . . . _^ gp-l ^ ^9-1. 
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Boundaries of the product domains can similarly be constructed. dS'^ is of course 
'zero'. Also, formally, we have Q, * dT^ ~ O for any fi. This follows from the requirement 
that O * 5T°, which is one of the boundary segments of * (the other being 50 *T^), 
agree with J7 itself. Thus, dT"~^ = OT^*" ~ nT"'^'^, that is n copies of T"—'^ joined at their 
boundaries. In the same language, we have d (S*^^^ * T'^^^'j ~ qS^~^ * T*?"^. 

Spherical domain T"~^, that can be viewed as constructed out of n copies of as T'^*", 
has n boundary segments of type T""^ that are perpendicular to each other, that is, have 
90° angles between their normals. It can be generalized to Tp~^, a similar but a 'distorted' 
domain, having angles cos~^ pij between the normals to their boundary segments i and 
j. A specific case is a regular Tp~^, denoted 7'(p)^^7 that has all those angles identical to 
say cos~^ p. As long as the n normals remain linearly independent, the n x n matrix p of 
Pij's completed with unity along the diagonal will be positive-definite. Writing p as AA^, 
Aik being the k^^ component of the i**^ normal in the a;-coordinatc system, and going to 
coordinates y = Ax, one obtains a simpler description of the Tp^^ domain in terms of that 
of T^~^, as one having the conical region M"q with the boundaries determined by = 
for the i*^ boundary segment. Inner products in the £C-system, such as x'^x' of two vectors 
X and x', can be computed in the y-system as p~^y' . 

5 Spectrum On The Sphere 

Many results are known in general about the spectrum of the Laplacian on compact domains. 
As already stated, all eigenvalues are real, non-negative, unbounded and discrete. The first 

eigenvalue has multiplicity one and the corresponding cigcnfunction can be taken to be 
positive within the domain. When the domain is the whole of the sphere S"""^, the resulting 
spectrum is well-known. In this case u is an integer taking values from zero to infinity. The 
first V, denoted z^i, is zero corresponding to a constant function on 5""^. The multiplicities 
of the eigenvalues will be revisited below. 

On the domain T"^-*^ under Dirichlct boundary condition, it is straightforward to show 
that ui = n. In fact, the simplest homogeneous function solving the Laplace equation 
in its cone M>q and vanishing on the boundaries is of degree n and is given simply by 
the product of the n-coordinates consistent with the first term in the expansion of sinh 
in equation (7). It is further clear that extending any spherical domain Jl"^^ by via 
a *-product would increase vi by one. If extending with instead, vi would remain the 
same. These observations are not trivial from the point of view on the sphere, but are a 
simple consequence of the property that z^i's are additive in a *-product. Additive property 
of z^i's in a *-product follows from the t~^'^~^ factor in the leading term (21) of the series 
solution and the multiplicative property of the heat kernel in a Cartesian product. 

To say more about the spectrum of the Laplacian on the sphere, let us next derive a 
spectral function, a 'generating function' for the eigenvalues and multiplicities in terms of 
f{x, x', t). Towards this end, let us set x' = x and r = 1 to obtain 

fix, X, 1) = r^-^e-^' ^^+^ (^') E(^-(^))'- (22) 

Note that a further operation of integrating over x, along with any f-independent weight, 
would integrate (h^a-i'f^))'^ to unity (its normalization) introducing the multiplicity rrii^. This 
procedure derives the following expression for the spectral function M(z) making use of the 
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Laplace transform of ij,, 



M{z) = V m^z^ = (1 - z^) z-"^ [ d^x e-^^^-^^^'-'fix, x, 1) 

-Inn 



(23) 



where < z < 1 and r is the length of x as given by (1). If the right side can be computed, 
this would provide us with both the eigenvalues and the multiplicities. 

The above function arose naturally from the solution of the heat equation on the cone. 
It differs from the usually studied trace of the heat kernel, Tre*^s, in that it is not the 
eigenvalues u{y + n — 2) of — V| that appear in the exponents, but rather z/'s themselves. Its 
derivation did not assume any specific character of the domain or the boundary conditions, 
except that D"- is conical intersecting S*"^^ into some domain ^V^~^. But its applicability 
depends on our knowledge of f{x,x,l). This is not expected to be the case in general. 
Below, let us first consider some special cases for which we do know /(a;, a;, 1). 

Consider again for the spherical domain the whole sphere S""^. In this case the associ- 
ated cone D'^ becomes the whole of M". Knowing f{x, x, 1) = (27r)~ 2 from (6), one readily 
obtains 



M{z) = (1 - z'^){l - zy = l + nz + J2 



k=2 



n + k 
k 



n + k — 3 
k-2 



(24) 



This gives the right eigenvalues and multiplicities on the whole sphere S*"^^. The two 
terms inside square brackets are the dimensions of the spaces of degree k and degree k — 2 
homogeneous polynomials in n variables, and the role of 1 — is hence to choose the 
difference for the dimension of the space of degree k harmonic homogeneous polynomials, 
that is those satisfying the Laplace equation in n-dimensions. 

Next consider the spherical domain T""^. In this case, under Dirichlet boundary con- 
dition, f{x, X, 1) is given by (7) that gives 



M{z) = z"(l - z^) 



2\l-n 



k=l 



n + k — 3 
k-1 



^n+2fe-2 



(25) 



For Neumann boundary condition, factor in front is absent. Extended domain Tp ^ for 



2, Tj^^j is a circular arc and, as 



p 7^ is not solvable similarly, except when n = 2. For n 

we have noted earlier, z/'s are multiples of ^ where = cos^^(— p) and are all of multiplicity 

— — 1 001 

one. Its Dirichlet spectral function is hence 2;'^(1 — z'fi)~ that becomes z {1 — z )~ for 
p = 0, corresponding to a quadrant arc in agreement with (25). 

An important property of M{z) that can be noted from its expression (23) is that, except 
for the factor 1 — z'^, it is nniltiplicativc in a *-product as a consequence of the multiplicative 
property of the heat kernel in a Cartesian product. Hence, M(z) on a spherical domain 

ns 

2 

1 



f2" ^ that is a *-product of the form $7"^ ^ * fig^ ^ is given by 



M{z) 



l-z- 



:Mi{z)M2{z), 



(26) 



where Mi{z) and M2{z) are the spectral functions on domains and ilg^"^ respectively. 

It follows from above that z^i's are additive in a *-product (as already noted). As an example 
of the above product rule, consider S^~^ * j"!-^ for which we obtain explicitly the Dirichlet 
spectral function 



M{z)=zi{l-z)-P{l-z^y-i. 



(27) 
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For instance on the half-sphere on S"'~^ that can be constructed as S'""^ * T^, wc get 
M{z) = z{l — zy~"' which incidentally has the same multiplicities as on T""^. Explicit 
results can also be obtained with one or more T^^^ 's included. Factor z'^ in front is absent 
for Neumann boundary condition. The above suggests the 'atomic' results: M{z) = 1 + z 
on S^, z on Dirichlet and 1 on Neumann T^. 

It is in general not straightforward to obtain explicit expressions for the eigenfunctions. 
However, the first eigenfunction hv^{f ) (of multiplicity one) of a *-product domain can be 
obtained from those of its factor domains via the multiplicative property of the following 
degree-i^i homogeneous function of a; in a Cartesian product (i^i being additive): 

KA^) = , r^'h^Ar)- (28) 

/r [y, + §) 

This follows from the leading term of the scries solution (21) and the multiplicative property 
of the heat kernel in a Cartesian product. On we have z^i = and hi^^{x) = 1, while 

on Dirichlet T""-*^ we have z/i = n and hu^{x) = 2^ IlILi suggesting the 'atomic' results: 
v\ = 0, h^^ (x) = 1 on 5° and ui = 1, hn-^ {x) = 2x on Dirichlet T^. If interested in exploring 
the hi,fj{f ) functions in general, we could rederive our spectral function without the angular 
integration to obtain 

/•oo 

M(f, f', z) = y\ m^{r, r')z'' = (1 - z'^)z-^ / dr r'^-'^e'^A^-'^ ^ f{rr, rf', 1), (29) 

Jo 

where m„{r, f') = J2fj hi,^{r)hi,^{r') . This provides us with a spectral function for the pro- 
jections on to the eigenspaces. As a function of zr with z considered as a radial coordinate, 
it can be identified as a kernel satisfying the Laplace equation on the cone for z < 1 
tending to 6{r — f') as 2; — > 1 under the considered boundary condition on In the case 

of M" having the whole sphere S""^ as the spherical domain, we get 

2 — 

where 6 is the angle between r and f', and jS*""^] is the size of the sphere S*""^ (surface 
area if 5^). This is the Poisson kernel of the n-dimensional unit ball at points zf and f 
that when expanded in powers of z gives rise to zonal harmonics as projections in terms 
of Gegenbauer (ultraspherical) polynomials. More involved expressions can be obtained in 
other explicit cases such as W x M?,q. 

One may also be interested in inverting (23) to obtain information about the heat kernel 
on the cone when the spectral function on a domain on the sphere is known. Given M{z) 
on a domain QP'"^, one can obtain fcin-i{t), where fnn-i{r^) = r"~^ fan-i d"'~^r f{x,x, 1), 

^ n 

as the Laplace inverse of fQn-i (s) = y^M{z),z = 1 + s — ^/s{s'+2), or obtain its 'trace' 

Jq drrfQn-i{r'^) as the Laplace inverse of ^/r2"-i(s). For instance, on the spherical do- 
main T""^, knowing M{z) from (25), fnn-i{r^) can be obtained as the Laplace inverse of 

2^ \ s(s+2) ^~^ ) ■ inverse is easily carried out for n = 2 to give 



1 1 , / ox .2 1 



dr fix, x,l) = ^- ^/o (r^) e-^ + -e-'^ . (31) 
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Here Qq is the quadrant arc and Iq is the modified Bessel function of order zero (this can 
also be obtained directly from the n = 2 series solution; alternately, knowing M{z), Laplace 
inverse can be viewed as summing up certain series of Bessel functions). More generally, 

one can obtain r"~^/(rr, rf', 1) inverting (29) as the Laplace inverse of j^M{f,r',z). 



6 Analytical Properties 

On continuing from the z < 1 region, M{z) exhibits a singularity at z = 1. At least for 
the various cases considered, the singularity is a pole of order n — 1 (the dimension of the 
sphere) so that we may write around z = 1 

^^(^' = (r^ + (r^ + -- <^^' 

Coefficients cq and ci can be determined, 

1 1 1(917"- H 

'^ = '^JS^\' ^1 = "2^° ~ 2 • ^^^^ 

It is convenient to write ci = —5(1 + 7)00 introducing 

Here, IS""^! and jS*""^! are the sizes of n — 1 and n — 2 dimensional spheres of unit radii 
respectively, lil"-"*^! is the size of the domain Jl""^ and |(9il"-"'^| is that of its boundary 
Sizes of and (90""^ are measured in units set by the n — 1 dimensional sphere 
S'^"^ of unit radius on which they reside. The sign of 7 is positive for Dirichlet and negative 
for Neumann boundary condition. Consistency with the product rule (26) implies that ^cq 
is multiplicative and 7 is additive in a *-product. 

The leading coefficient cq can be determined by letting z 1 in the expression for M{z). 
Note that the exponential inside the integral would no longer provide the suppression as 
r ^ 00. As r — >■ 00, f{x,x, 1) away from the boundary tends to a constant (27r)~2 (see 
equation (6) at a; = a;', r = 1). The integral is thus dominated by regions near r = 00 where 
the angular integral contributes |ri"~^|. This gives, ase = l — 2:^0, 

M(l-e)~2e/ dr r"-V^^'^'J J = 2-!^^ J ^. (35) 

^ ' Jo (27r)i 27ri e""^ ^ ^ 

The factors in front can be identified as twice the inverse size of the sphere S'^~^. 

The next coefficient ci can be determined by the method of images. To start with, note 
that the contribution to M{z) coming from the source alone, 

Cq 1 + Z ^ Cq 1 Cq 

2(1- z)"-i ~ (1 - z)"-i 2(1- z)^-^ ' ^ ^ 

makes an order n — 1 contribution as well. In the method of images, the source placed 
within the domain induces images across the boundary that, under Dirichlet condition, 
cancel out the source effect on the boundary to ensure zero boundary condition. Since 
/(£C, £C, 1) is evaluated at the source location itself, as a; is varied, the source moves and 



9 



the images follow the source. As r — t- oo many of the images will recede away from the 
source. The leading contribution comes from the image brought closest to the source by 
taking the source close to the boundary. Its contribution is ~ — (27r)^ 2 e~^^-L. Here y±_ 
is the perpendicular distance of the source to the boundary so that the image to source 
distance is 2yx ■ The image contribution ase=l — z^Ois 

-T^TTT / drr^ 2- / dy^e = -- ' ' ■ (37) 

(27r)2 7o ian± ^ 27r~ ^ 

The factors in front can be identified as half the inverse size of the sphere A negative 

sign is chosen to satisfy Dirichlet condition on the boundary. For Neumann boundary 

condition, the sign will be positive. 

Expansion (32) is a result of an expansion of /(a;, a;, 1) in r^^ in the expression (23) for 
M{z). Since f(^x,x,T) is function of the combination an expansion of f{x,x,l) in 

is in fact an expansion of f{x,x,T) in at r = 1. This is the well-known expansion 
of the heat kernel (see for instance Vassilevich [2003]), in our case on the cone D"'. Because 
the higher order terms of this expansion bring in more powers of r into the denominator 
inside the integral in (23), as such it can only be used upto coefficient Cn-i- If the remainder 
falls off faster than r"" as r — >■ 00, its integral will be finite at 2; = 1 because of the r — >■ 
behavior of f{x,x,l) evident from (22). Also note here that the heat kernel expansion 
being an expansion in does not see any terms of the type e^'' for instance. That such 
terms are present can be seen by taking the example of the n = 2 domain for which we 
know J^i dr f{x, x, 1) from (31). The first two terms on the right hand side of (31) give rise 
to the heat kernel expansion while the last term, not visible to the heat kernel asymptotics, 
is required for the r — t- behavior. 

Expansion (32) can also be obtained from the heat kernel expansion on O""^ on the 
sphere itself. This can be done using the identity^ 



,x se'^" dt ,v2 e^* f°° dt ^ Ps^^ ^ .2 

M(e-^) = ^^ —e-^^-4iTie^^s= — e-*-TrTre«^s, (38) 

2V7r Jq t2 V7I" Jo Vt 

where I = ^(n — 2) and Tre*^s is trace of the heat kernel on the sphere. The identity follows 
from the representation of the trace as ^^miye~^'^^'^^'^^K Here it is convenient to expand 
M (e~*) around s = 0, 

M(e-)c. + + + (39) 

The heat kernel expansion on the sphere can be expressed as 

Tre*^s c± {ATrt)-"^ (qq + aiVt + a2t + ■ ■ (40) 

Coefficients of this expansion are known to certain order (see for instance Vassilevich [2003]), 
the first three of which are 



ao = |0"-i|, ai = T^|50«-i|, as = J / R+lf K. (41) 



^Analogous relation can be written down for the pointwise object M{r,r' ,z). Inverse relations can be 
obtained by expressing them as Laplace transforms, giving rise to identities for the heat kernel such as the 
one involving the Jacobi 0- function on . 
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Here R = (n — l)(n — 2) is the scalar curvature of S"'~^ and hence that of $7""^, and K is 
the trace of the extrinsic curvature of the boundary dil'^~^ relative to The sign of ai 

is negative for Dirichlet and positive for Neumann boundary condition. Use of the above 
expansions in the identity (38) gives us the 6-coefHcients, 

where £ = ^{n — 2) and, cq and 7 are given by (33) and (34). The expression (41) for 02 is 
applicable on domains with smooth boundaries (as well as on those without boundaries). On 
domains that have corners on their boundaries, it is applicable with additional contributions 
coming from the corner regions. For instance, on a domain on S'^, each vertex of angle ip on 

the boundary would contribute g (^■^ ~ to 02 (making use of a result attributed to Kac 

[1966] that can be inferred in the present framework from the trace of the heat kernel on 
the n = 2 cone following the discussion at the end of section 5). On a domain on a higher 
dimensional sphere, this result is applicable when integrated along the corner regions of 
dimension n — 3. For instance, on T"~^, one finds that the corner contribution to 02 is 
^n(n — l)(n — 2) |T"~^| (contribution from K being zero) consistent with the expansion of 
the explicit result (25). 

The series expansion of the kind (39) is useful in estimating the growth of the spectrum 
at large eigenvalues. This is done with the help of a counting function 

W(z/) = ^m,.W, (43) 
I/' 

where l,^'<,y is the step-function. W{u) counts the eigenvalues including multiplicity up to 
u. Its Laplace transform is 

W{s)= r dvW{-u)e-"' = -M{e''). (44) 
Jo s 

As we have noted, M{e~^) is expected to have a pole of order n — 1 at s = giving rise to 
an expansion of the kind (39). Here the singularity should arise from the large v behavior 
of W{y). One finds for W{y) the large v expansion 

f, ,,n— 1 I. ,,n—2 u ,,n— 3 

m-)~7^ + 7^ + 7^ + ---- (45) 

(n — 1)! (n — 2)! [n — 6)1 

Expressed in terms of the eigenvalues A = z/(z/ + n — 2) ~ i/^ of the Laplacian on Q""-^, this 
is consistent with the Weyl scaling law (true for more general domains). 

As a Dirichlet series in s, one expects M(e~*) defined on the positive real s-axis to 
be analytic on the half-plane Re(s) > 0. Its behavior for Re(s) < is less clear. Result 
(23) indicates naively a relation M = — z"~^M(z). However, this is not expected 

to hold as an approach to from z along the real axis encounters the singularity at 
z = 1. For the cases considered earlier, one finds instead M {z'^) = {-1^-^ z''-^-'' M(z) 
as well as Md [z'^) = {—1)'^~^z"'~'^Mn{z) where subscripts refer to Dirichlet and Neumann 
boundary conditions. Being consistent with the product formula (26), these will also hold 
for domains factorizable into such cases. They are however restrictive to hold in general, 
but when one does, M(e~^) can be expected to be analytic on the half-plane Re(s) < 
(with singularities along the imaginary s-axis). 
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7 A Scaling Procedure 



It is a result that the eigenvalues of the Laplacian do not increase as the domain is en- 
larged. Having dimensions of inverse coordinate squared, eigenvalues can be expected to 
scale accordingly, though in general approximately, suggesting that we look for a scaling 
procedure to estimate the eigenvalues on a spherical domain. However, applying scaling 
to the eigenvalues itself, as is usually done, turns out to be not satisfactory. Let us hence 
look for a spectral function M[z) on a target domain Jl""^ of the form (under Dirichlet 
boundary condition) 

M{z) = z"Mo(z^), (46) 

where Mq{z) is the known spectral function on a reference domain f^o~^- This implies that, 
given the eigenvalues Aofe = t'ofe(^ofc + n — 2),k = 1, 2, • • • of the Laplacian on fig"^, the 
eigenvalues = Vki^k + n — 2) on can be estimated according to 

Vk = a + l3vok, k = l,2,---. (47) 

Parameters a and /3 can be determined by expanding M{z) and Mq{z) into their series (32) 
at z = 1 and matching the first two coefficients (33) for the two domains, 

1 

n-l 

(48) 

where 7 and 70 for and are as given by (34). This estimation procedure can 

also be expressed as a scaling of the combination z/ + ^ (n — 2 — 7) . Note that this does 
not change multiplicities. If and ^q~^ are closely related and the eigenvalues are 

well separated, this may be a reasonable assumption to make; at least for the first few 
eigenvalues. Eigenfunctions will of course be different. 

The linear scaling estimate (47) with its two parameters is able to match the first two 
coefficients of the series expansion (32). It is helpful to have it extended to involve three 
parameters to be able to match the third coefficient of the series as well. A convenient 
approach to deriving such an extension without affecting multiplicities is to require the 
counting function W{y) to agree for the two comparison domains. Keeping the first two 
terms of the expansion (45), W(y) can be approximated for large v SlS ^ {n-i)\ + 

with p = To this order, one then obtains the scaling procedure (47) as a scaling of 
the linear combination v + p. To include the next order term, let us look for a quadratic 
combination {y + p)"^ + q for some p and q such that for large u 

Expanding this, comparing with (45) and using the 6-coefficients from (42), one obtains 

P = ^-l^ Q = -^'-li^-^h' + 0zry (50) 

where ^ = ^(77, — 2), 7 is as defined in (34) and 02 is given in (41). Thus, given p and q as 
above, {i^ + p)^ + q for a domain Jl""^ can be estimated as times the same combination 
for a reference domain Scaling factor ^ is the same as before given by (48). 
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The quadratic scaUng procedure thus defined is quite general apphcable to domains with 
smooth boundaries under Dirichlet boundary condition. It can also be used on domains 
without boundaries in which case it becomes a linear procedure applied to the eigenvalues 
of the Laplacian itself. On domains that have corners on their boundaries, the procedure is 
applicable with additional contributions to 02 coming from the corner regions as discussed in 
section 6. Under Neumann boundary condition, the scaling procedure needs to be modified 
so as to preserve ui = 0. The equivalent of the linear procedure in this case is to scale 
^{ly + 2p) by The equivalent of the quadratic procedure turns out to be to scale {v + 
pY + ^qu — by Both p and q as well as /3 are given by the same expressions as those 
obtained above (note however that 7 is now of the opposite sign). 



8 Numerical Comparisons 

The following numerical comparisons are for domains on the two-sphere of unit radius 
under Dirichlet boundary conditions, carried out with the linear scaling procedure (47). 
For clarity, area is denoted as A and the perimeter \dVl'^\ as L so that for the scaling 

parameters (48), we have 7 = ^, 70 = ^, /? = ^J^, a = ^[7 - 1 - /3(7o - 1)]. 

For the domain T^^-^, the domain size expression is given by (60). As p ^ 1, Tj^^^ 

tends to cover half the sphere. Since M{z) = z{l — z)^"" of the half-sphere is exactly 
related by scaling to Mq[z) = z^{l — z^)^^" of the p = domain T"^^, scaling estimates 
can be expected to be reasonable for p in-between. Domain T^^^ on the two-sphere is a 
spherical equilateral triangle of vertex angle cos^^(— p). It has A = 3cos^^(— p) — vr and 
L = 3 cos~^ (~T+p) • ^''^ case, the reference domain for scaling estimation can be chosen 

to be the octant triangle having Mo{z) = z^{l - z'^)~'^, Aq = |,Lo = ^,70 = t'oi = 3. 
For p = ^, the spherical equilateral triangle is a tetrahedral triangle having A = vr, L = 
Scos"-*^ ("5)- R-atzkin and Treibergs [2009], in the context of a capture problem, present 
a theoretical and numerical result 5.159 for the first eigenvalue Ai = J^i(i^i -|- 1) of the 
Laplacian on the tetrahedral triangle. Scaling estimate gives vi = 1.826 and Ai = 5.162 in 
excellent agreement with their result. 

A spherical cap is a disk-like domain on the two-sphere. If its radius relative to its center 
in angles is 9, it has A = 2it{1 — cos 9) and L = 27rsin^. In this case the reference domain 
can be chosen to be the half-sphere that has Mq{z) = z{l — 2;)"^, 70 = 1^01 = 1 so that 

The usual scaling procedure applied to the eigenvalues of the Laplacian itself is based on 
just the size of the domain, and hence is not able to differentiate the effects of the boundary. 
Ratzkin and Treibergs [2009] present a theoretical result Ai = 4.936 for the first eigenvalue 
on a spherical cap (^ = |) having the same area as the tetrahedral triangle. Scaling with 
(51) gives Ai = 4.949 in excellent agreement. 

A sector of the spherical cap making an angle at its center has A = Lp{l — cos 9) and 
L = 99 sin -|- 29. Choosing the reference domain to be such a sector on the half-sphere (a 

product domain T^_(.Qg^^ * T^) that has Mo{z) = z^~^^{l — z'^)~^ ^1 — zv^ ,70 = foi = 
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1 + f > we get 



^1 = o cot - + - ^ . , - 1 + (52) 



2 v^sin^f \/2(^sin| / \/2 



sin ^ 



Ratzkin and Treibergs [2009] present a theoretical result Ai = 5.0046 for the case (p = ^ 



2n 



and 9 = cos ^ whereas the scaling procedure gives Ai = 5.1046. 

As a domain on the sphere is shrunk retaining its shape, it tends to approximate a flat 
domain in the limit, allowing for a comparison to the available solutions on flat domains. For 
instance, as the spherical cap has its radius 6* ~ 0, its z^i ^ (1 + ^2)5"^ = 2.4142(5"^ 
that compares well with the flat disk solution ^/X^ = jo,iS''^ = 2.4048(5"^ (jo,i being the 
first zero of the Bessel function Jq). The second one 1/2 ~ (1 + 2^2)5"^ = 3.8284(5-^ also 
compares well with = ji.i^^^ = 3.8317(5^^. The next one ^3 ~ 5.2426 is close to 
= 5.1356. As expected, higher ones start showing up significant differences. 

Complete solution on the equilateral triangle on the plane was obtained by Lame [1833]. 
Comparing the octant triangle on the sphere, one finds for the equilateral triangle of side 
length S on the plane vi ~ ^2\/3 + 2y^^^ = 7.2736"^ that compares well with Lame's 

result = ^6-^ = 7.255S-\ The second one 1/2 ~ (2^3 + 4^^^) S'^ = 11.083(5"^ 

also compares well with -v/A2 = ^^^^S~^ = 11.0825"^. The next one 1/3 ~ 14.892 is close to 
s/Xs = 14.510. Here too, higher ones start showing up significant diff'erences. 

More generally, one can use the flat domain solution to estimate the first few eigenvalues 
on a similar domain on the sphere. Given A and L for a domain on the sphere and Aq 
and Lq for a similar domain Jig on the plane, one finds for r^^, = 1, 2, • • • on O^, 



1 fL Lq \ ApXok 



2 VAA^ J \l A 

where Aofe, k = 1, 2, • • • are the eigenvalues of the Laplacian on CIq. This may be viewed as 
providing a curvature correction to the flat space eigenvalues. Note that the length scale 
on Qq cancels out, and that A, L and 1^1 are in units set by the unit sphere. 

The above comparisons are done using the linear scaling procedure. Better results can 
be obtained with the quadratic scaling procedure discussed before. The quadratic procedure 
can be used on domains with smooth boundaries, for instance on a spherical cap on S*""^ 
(K = (n — 2) cot^, where 9 is its radius in angles), but it turns out to be identical to the 
linear one in the case of the spherical cap on S^. As noted earlier, the quadratic procedure 
is also applicable on domains that have corners on their boundaries. For our numerical 
comparisons, this offers the improvements: Ai = 5.1625 — t- 5.1606 for the tctrahedral 
triangle, Ai = 5.1046 — 5.0187 for the sector of a spherical cap and ^/Xl = 7.27345"^ — >■ 
7.26135"^ for the equilateral triangle on the plane. 

Viewed as an extension of the scaling of the Laplacian eigenvalues on flat domains to 
include curvature effects, the discrepancies in the scaling estimates could become significant 
as domains get too large relative to say the half-sphere. Also, the scaling procedure based 
on just two or three parameters is not expected to yield good results for all the eigenvalues, 
but its potential to do so for the first few is intriguing, especially because it is based on 
the first few coefficients of the scries that governs the growth of the spectrum at large 
eigenvalues. It will be interesting to study the applicability of a similar scaling procedure 
to the spectrum of other differential operators, or to more general domains extended to a 
cone or by taking (38) as defining M{z). 
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A Spherical Domain Sizes 

Size of a spherical domain denoted can be computed as the t = 1 value of 

Z{t) = 2 [ (rxS{r'^-t). (54) 
Jd" 

A factor 2 is introduced since S{r'^ - 1) = l6{r - 1). The integral over X reduces to an 
integral over the angular variables giving rise to the size of the spherical domain at t = 1. 
Laplace transform of Z{t), denoted Z{s), is 

Z{s) = i / d^xe"'"'. (55) 

Conical nature of allowed us to scale x by s~2. Inverting the transform, we get 

|0"-i| = |5"-i|4r / d^xe-'\ |5-i| = -^. (56) 

This gives us the size-fraction taken as a fraction of the size of the sphere of the same 
dimension. As can be noted from above, size-fractions arc multiplicative in a *-product. 
For the spherical domain T"^^, = IR"q and the above expression gives 2^" for its size- 
fraction consistent with its description given earlier. This also follows from the multiplicative 
property of size-fractions given that, in this setup, we have l^^l = 2 and |r°| = 1 so that 
the size- fraction of can be taken to be ^. Similarly, the size- fraction of S'*'"^ *r«-i turns 
out to be 2^"? as expected. 

As noted in section 4, spherical domain T""^ can be generalized to T'p^^ having elements 
of a matrix p for the cosine of angles between the normals to their n boundary segments. 
Writing p as AA^ and going from coordinates a; to y = Ax, we obtain from (56) the domain 
size expression 

I on— 1 1 f 

|t;-i| = 4^-^/ <Pye-y''p-'y. (57) 

7r2 ^/detp J«^^^ 

Region of ^/-integration is M"q since the boundaries are now determined by t/j = for the fi^ 
boundary segment for each z = 1, • • • ,ra. |0T"~^| can be computed from the above formula 
for each of the boundary segments, with p restricted to one dimension less. 

An example of Tp~^ is a regular domain TJ^^^^ with a single parameter p having angle 

cos~^ p between the normals to their boundary segments. It has the p-matrix 

P., = (i-P)5.,+P, Pi.— r^^^.- (i_^)(i + (,_i)^) - (58) 

This has determinant detp = (1 — p)"~^(l + {n — l)p). Matrix A can be chosen to be 

Ay = a5ij + 6, a = ^l- p, b=^ (^^/l + {n-l)p - ^1 - p) ■ (59) 
For this domain, the domain size expression (57) simplifies to (p > 0) 



rpn — l 



1 



= |r"-i| — / due-"" 



2 



erfc 



VP 



u 



(60) 
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where IT 



•n—l I 



2 "' 5'" 1 and erfc is the complementary error function. The same expres- 



sion upon setting n 



n 



1 and p ^ gives ^ 



for any n giving 



nn— 1 
■(1/2) 



1 5" ^1 corresponding to a domain on 5" ^ analogous to a 



It can be evaluated for p 

:in— 1 



tetrahedral triangle on S . For general p, on obtains the recursive differential equation 



n[n — 1) 



fn-2 



l + 2p 



Up) 



pn— 1 



^n-l\ 



(61) 



= cos ''^(— p) as expected 
Tj)| = 3cos-i(-p)-7r 

as expected for a spherical equilateral triangle of vertex angle cos~^(— p) (and also having 



with fo{p) = flip) = /n(0) = 1. Solving forn = 2 gives 
for a circular arc T^^^ of angle cos~^(— p), and for n = 3 gives 



Scos-i (-TT^))- For small p, /„(/9) ^ l + ^n{n-l)p+^n{n-l){n-2){n-3)p^. 



Since fn{p) is an increasing function of p, we have 



so that T^"^ ^ tends to cover half the sphere. 



t; 



■n— 1 
(P) 



> \T 



'n— 1 1 



. As p ^ 1, Up) 2 



n-l 
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